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_ We present a study performed for the GLIG CDR on the measurement of chargino 

and neutrahno production at ^/s = 3 TeV. Fully hadronic final states with four jets 
(N- and missing transverse energy were considered. Results obtained using full detector 

simulation for the masses and for the production cross sections of the changino and the 
lightest and next-to-lightest neutralinos are discussed. 

Qj" 1 Introduction 

oT 

i-S^ ■ For the CLIC Conceptual Design Report (CDR) [I], several physics processes were selected to 

benchmark the performance of general purpose detectors at a centre-of-mass energy of 3 TeV. 
The study presented in this document assumes a SUSY model defined by the mSUGRA 
parameters mi/2 = 800 GeV, Aq = 0, mo = 966 GeV, tan /3 = 51 and /z > 0. In this model, 
the lightest chargino, xt ^ has a mass of 643.2 GeV, while the masses of the lightest and 
next-to-lightest neutralinos, Xi £^nd ^^'^ given by 340.3 GeV and 643.1 GeV, respectively. 
The lightest neutral Higgs boson, has a mass of 118.5 GeV. 

The pair production of charginos and neutralinos was investigated: 

e+e- ^ xf xr ^ W+x"iW-x"i and (1) 

e+e-^x'2x'2^h'iZ')x>'iZ')x'i, (2) 

where BR{xf VV^Xi) = 100%, BR{x^ h^xl) = 90.6% and BR{x^ Z°x°^) = 9.4%. 
Hadronic decays of the W^, hP and Z° bosons were considered and hence the investigated 
final state is given by four quarks and missing transverse energy. The reconstruction of 
chargino and neutralino pair production allows to benchmark the reconstruction of hadron- 
ically decaying gauge bosons in multi-liadron final states. 

More details on the study discussed in the following are given in [2]. 



2 Monte Carlo production 

The physics events used for the study presented here were generated using the WHIZARD 
1.95 [3] program. Initial and final state radiation (ISR and FSR) were enabled during the 
event generation. The luminosity spectrum expected at CLIC was used during the event 
generation. The hadronisation of final state partons was simulated using PYTHIA [4]. The 
generated events were accordingly passed through the detector simulation program SLIC 
which is based on the Geant4 [5] package. The CLIC_SiD_CDR [6] detector geometry model 
was used. 
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Events were overlayed with pileup from 77 — > hadrons interactions corresponding to 60 
bunch crossings [7] . The reconstruction chain included an improved version of the Pando- 
raPFA [8] algorithm to reconstruct particle flow objects. 



Type 


Process 


Cross section [fb] 


Luminosity [ab ^] 


Referenced with 


Signal 


xtxl 

X'2X°2 


10.6 
3.3 


13.4 
23.8 


Chargino 
Neutralino 


Background 


xtX2 
xtx2 
XiXi^v 

X2X2'^'^ 


10.5 
0.8 
1.4 
1.2 


1.8 
8.9 
21.9 
13.3 


SUSY 




qqqqz^F 
qq/i°z^F 


95.4 
3.1 
0.6 


4.5 
6.2 
22.8 


SM 



Table 1: Cross sections and integrated luminosities of the available Monte Carlo samples for 
chargino and neutralino pair production and for SUSY and Standard Model backgrounds. 
The charge conjugated modes are implied throughout this document. 



An overview of all produced Monte Carlo (MC) samples is given in Tab. [TJ Dedicated 
samples for the considered signals corresponding to large luminosities are available. Ad- 
ditionally, backgrounds from SUSY and Standard Model (SM) processes were used in the 
study presented in this note. 

3 Event reconstruction 

The steps to reconstruct events with four jets from particle flow objects (PFOs) are described 
in this section. The presence of pileup from the process 77 — !■ hadrons increases the number 
of reconstructed PFOs in typical signal events by a factor 10 and the total visible momentum 
by a factor four. On the other hand, the background particles are emitted mostly in the 
forward direction. A large fraction of the background was rejected using combined timing 
and momentum cuts. 

Only events containing at least four reconstructed PFOs with px > 250 McV were 
considered further. Events with at least one identified electron or muon with px > 20 GeV 
were rejected. 

Jets were reconstructed from PFOs using the kt algorithm [9] as implemented in Fast- 
Jet [TU] in its exclusive mode with R = 0.7 and using the E recombination scheme. The 
clustering was stopped when four jets were found. To reject leptonic decays of W^, or 
Higgs bosons further, all jets were required to contain more than one PFO. 

Bosons candidates were formed from jet pairs minimising: 

(Mjj- 1 - Mw±.hof + (Mjj,2 - Mw±,ho)^, (3) 

where Mjj^i and Mjj^2 are the masses of the two reconstructed jet pairs and M^r±^f^o was 
set to the world average of the boson mass to reconstruct xf and to the assumed Higgs 
boson mass to reconstruct 
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The reconstruction of bosons in XiXi events is illustrated in Fig. [T] The distri- 
butions obtained with and without the overlay of 77 hadrons are compared. A good 
reconstruction of bosons was achieved if combined timing and momentum cuts were 
applied to select the PFOs used as input to the jet reconstruction. The right plot in Fig. [1] 
shows the reconstructed masses for the two selected signal data sets: chargino pairs decaying 
into W'^W~ and neutralino pairs decaying into either hPhP or hPZ'^ . Since less than 1% of 
the neutralino pairs decay to the Z^Z^ final state, this contribution is not shown in the fig- 
ure. The horizontal band for Mjj^2 ~ M^^a and Mjj^i < My^o is caused by X2X2 ~^ f^-'^Xi^'^Xi 
events where one of the h'^ bosons is only partially reconstructed. No similar vertical band 
is visible due to the way the jets are ordered in the analysis. 
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Figure 1: Reconstructed mass of W candidates without and with pileup from 77 — > hadrons 
(left). For events with pileup, the distribution for all PFOs is compared to that obtained 
using combined timing and momentum cuts ( "Tight selected PFOs" ) . 2D mass plot of the 
signal final states for different decays of chargino and neutralino pairs (right). The peaks 
corresponding to the individual chargino and neutralino decays are indicated. The event 
samples were scaled to have similar number of events for each channel. 



4 Event selection 

The selection of xtxi ^-i^d X2X2 P^i'" production events was performed in two steps. First, 
a cut-based preselection was applied. The remaining backgrounds were suppressed further 
using boosted decision trees in a second step. These two steps arc described in the following 
two subsections. 

4.1 Preselection cuts 

To restrict the training of the boosted decision trees to the region where the signal purities 
are high, the following preselection cuts were applied: 

• 40 < Afjj,i < 160 GeV and 40 < A/jj,2 < 160 GeV 
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• I COS 6'™'^^ I < 0.95, where 6*™'^^ is the polar angle of the missing momentum 

• Angle between the or Higgs candidates larger than 1 radian 

• I cos6'jj'^| < 0.95 and | cos0jj'^| < 0.95, where 9^''-^ and 6'^''^ are the polar angles of the 
two jet pairs 



4.2 Event selection using boosted decision trees 

To distinguish between signal and background events further, boosted decision trees as 
implemented in TMVA [TT] were used. For training purposes, 20% of the available events 
for each process were used. These events were not considered in the analysis to measure 
masses or cross sections. 

The boosted decision trees were trained using 15 variables describing the event topology 
and describing kinematic quantities of the reconstructed or Higgs candidates. The 
efficiencies of the entire selection chain consisting of the preselection and of the BDT cut for 
reconstructed chargino and ncutralino signal events were 25% and 33%, respectively. The 
signal purities in the selected samples were 57% for the chargino and 55% for the neutralino. 




Figure 2: Reconstructed energy (left) and mass (right) for the chargino signal and for 
the SM and SUSY backgrounds. The histograms are stacked on top of each other. All 
distributions are scaled to an integrated luminosity of 2 ab~^. 



The energy and mass distributions of reconstructed Chargino and Neutralino candidates 
are shown in Figs. [2] and [H 



5 Signal extraction 

Two complementary methods were used to extract masses and cross sections from the energy 
distributions of the reconstructed and selected and Higgs candidates (see Figs. [2] and [3]). 
Both approaches are described in the following subsections. 
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Figure 3: Reconstructed Higgs energy (left) and mass (right) for the neutrahno signal and 
for the SM and SUSY backgrounds. The histograms are stacked on top of each other. All 
distributions are scaled to an integrated luminosity of 2 ab^^. 



Parameter 1 


Uncertainty 


Parameter 2 


Uncertainty 




6.3 GeV 




2.2% 


Mm 


3.0 GeV 




1.8% 


Mm 


7.3 GeV 


<^ixm) 


2.9% 



Table 2: Uncertainties of the chargino and neutralino masses and pair production cross 
sections obtained from two parameter template fits. An integrated luminosity of 2 ab~^ is 
assumed. 



5.1 Template fit 

The pair production cross sections and masses of the xf and ^2 particles were determined 
using a template method where chargino and neutralino signal Monte Carlo samples for 
different mass hypotheses were generated. The M{xi) mass was also measured since the 
energy distribution of bosons from xf decays is sensitive to this observable. Two- 
dimensional fits were performed simultaneously to the mass and production cross section 
for a given particle to account for the correlation between both quantities. The statistical 
uncertainties of the extracted masses and cross sections, determined using toy Monte Carlos, 
are shown in Table [2] All measured values are in agreement with the input values used in 
the Monte Carlo generation. 

5.2 Least squares fitting 

Linear least squares fits of gaugino masses and cross sections were used to check the results of 
the template fitting technique and to fit for three or more parameters simultaneously. Each 
W^/Z°/h° reconstructed energy histogram bin was expanded linearly about the nominal 
masses and cross sections. The slopes were obtained by convoluting a map of true-to- 
reconstructed bin contents with the true energy distributions at different gaugino masses. 
No fits were actually performed; instead the bin statistical errors were calculated and then 
propagated to the fit parameter errors using standard formulae for linear least squares fits. 
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Parameter 


Error 


M(xf) 


7.3 GeV 


M(x?) 


2.9 GeV 




9.8 GeV 


<^{xtxi ) 


2.4% 


<^{X°2X°2) 


3.2% 



Table 3: Statistical mass errors and relative cross section errors for a five parameter least 
squares fit which includes the constraint that the Xi mass be within 3 GeV of the value 
measured from slepton analyses. 



The least squares results for two parameter fits are in reasonable agreement with the 
template fits described above. The two parameter fits assume that the other SUSY param- 
eters have been measured with arbitrary accuracy. The Xi mass will be measured with an 
accuracy of AM(xi) = 3 GeV at CLIC by combining the results from the slepton anal- 
yses [H]. A term constraining the Xi mass to be within 3 GeV of the best estimate was 
added to the least squares fit. One can combine the data from the xt and X2 analyses and 
perform a five parameter least squares fit of M{xf), M{x2), M{xi),cr{xiXi)T'^{X2X2)- The 
results for this fit are shown in Tab. El 



6 Summary 

The signals from xf and X2 pair production were extracted from fully hadronic final states 
with four jets and missing transverse energy. The study was performed using full simulation 
and considering pileup from 77 — > hadrons. Two different signal extraction procedures are 
in reasonable agreement. The chargino and neutralino pair production cross sections were 
extracted with a statistical precision of 2 — 3% while the masses of the xt^ Xi and X2 
particles were determined with typical statistical accuracies of about 1 — 1.5%. 
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